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We consider the Susceptible-Infected-Recovered (SIR) epidemic model on a Euclidean network in 
one dimension in which nodes at a distance I are connected with probability P{1) oc l~^ in addition 
to nearest neighbors. The topology of the network changes as 5 is varied and its effect on the SIR 
model is studied. R{t), the recovered fraction of population up to time t, and r, the total duration 
of the epidemic are calculated for different values of the infection probability q and S. A threshold 
behavior is observed for all 5 up to 5 « 2.0; above the threshold value q — Qc, the saturation value 
Raat attains a finite value. Both Raat and r show scaling behavior in a finite system of size A*'; 
Rsat ~ iV"'^''"ffi[(g - qc)N^^''] and r ~ N''^''g2[{q ~ qc)N^'"]. qc is constant for < 5 < 1 and 
increases with S for 1 < 5 < 2. Mean field behavior is seen up to (5 « 1.3; weak dependence on 5 is 
observed beyond this value of S. The distribution of the outbreak sizes is also estimated and found 
to be unimodal for q < q^ and bimodal for q > qc- The results are compared to static percolation 
phenomena and also to mean field results for finite systems. Discussions on the properties of the 
Euclidean network are made in the light of the present results. 

PACS numbers: 64.60. Cn, 64.60. aq, 64.60. F-, 64.60.ah 



I. INTRODUCTION 

Studies of familiar static and dynamic phenomena on 
complex networks have led to some surprises [IHSl in re- 
cent times. For example, the Ising model which in one 
dimension does not have a phase transition, showed not 
only the existence of a phase transition but also that it 
occurs with mean field criticality 0, on the Watts- 
Strogatz (WS) [1| type network. On the scale-free net- 
work also, it showed a behavior not encountered usually 
on regular lattices; the transition temperature showed a 
logarithmic increase with the system size Q • Among the 
well studied dynamical phenomena on networks IS* are 
opinion dynamics models ^10] , disease and damage prop- 
agation |ll| - [l4j . synchronization of coupled oscillations 
[2, EBl , zero temperature quench of Ising and other spin 
models [T6l-[20|. etc. Just like the static results on net- 
works, there have been special features of the dynamical 
phenomena as well. On complex networks, zero temper- 
ature quenching of Ising model shows freezing even in a 
one dimensional network with additional random links. 
Another example is the voter model: its dynamics can 
be conceived in different manners on networks leading to 
the conservation of different physical quantities. How- 
ever, on regular lattices these different dynamical rules 
are equivalent jlQI]. 

A standard model for epidemic spreading is the 
Susceptible-Infected-Recovered (SIR) model where indi- 
viduals can be in three possible states; susceptible: who 
are liable to be infected, recovered: those who contracted 
the disease but are now recovered and immunized, in- 
fected: people who are suffering from the disease and 
can infect others. A very well known result for the SIR 
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model is that it is like a dynamical percolation prob- 
lem and its critical behavior coincides with that of static 
isotropic percolation on regular lattices [2lj. On the 
other hand, another model for disease spreading, namely 
the Susceptible-Infected-Susceptible (SIS) is identical to 
a directed percolation problem. Percolation phenomena 
has also been studied extensively on complex networks 
(earliest works appear in [22| - [2^ and comprehensive re- 
views are available in ^'^'^ possible to compare 
the two phenomena of epidemic spreading and percola- 
tion on complex networks as well. 

Although extensive research work on epidemic spread- 
ing on scale free and WS networks have been made, it 
may be noted that social connections are neither scale 
free nor like those considered in the WS network. In fact, 
many social networks have a spatial dependence in the 
connection probability of the nodes [l^. We thus con- 
sider SIR on a spatial model in which random long range 
links between nodes at a Euclidean distance I are added 
with probability P[l) oc . Nearest neighbor links are 
always present in this model. In [26l , [27| , a similar model 
in two dimensions was considered, but the probability 
P{1) was essentially dictated by the heterogeneity of the 
degree distribution, and only short range links followed 
the power law distribution while long range links were 
added randomly. 

Static properties of the network considered in this pa- 
per are quite well studied [H, [2^ . The network behaves 
as a small world network for 5 < 1 and as a regular one 
dimensional lattice for 5 > 2. Some ambiguities remain 
regarding its behavior when 1 < 5 < 2. While some 
results suggest the network is still a small world here, 
other works claim that it has a finite lattice-like behavior. 
Dynamic phenomena like zero temperature coarsening of 
Ising model and searching have also been studied on this 
network 

The aim is to locate the infection threshold values as 
a function of 5 and find out the critical behavior for SIR 
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on the Euclidean network. This may also help in under- 
standing the nature of the network in the controversial 
region 1 < 5 < 2 by considering a dynamical process. We 
have also studied the way the recovered fraction grows in 
time as well as the total duration of the entire epidemic 
spreading process. 

In section II, we describe the model and method briefly 
and follow it up with the results in section III. In section 
IV, we discuss a simple model (mean field type) to make 
some comparisons. In section V, the results are discussed 
and compared to earlier studies made on the same net- 
work. In the last section concluding remarks are made. 



II. MODEL AND METHODS 

The network is generated in the following way: in a 
system of N nodes placed on the sites of a one dimen- 
sional lattice, all nearest neighbors are first connected. 
Additional long range links (one per node on an average) 
are established at a distance Z > 1. This is done in the fol- 
lowing way: two nodes are selected randomly; if they are 
not connected already, a connection is established with 
a probability P{1) cx {I is the distance separating the 
nodes). The process is stopped as soon as N/2 links have 
been formed this way. The average degree of the nodes 
is therefore three. 

In a homogeneous system, the SIR model can be de- 
scribed in terms of the densities of susceptible, infected 
and recovered nodes, S{t), I{t) and i?(i), respectively, 
as function of time. These three densities are related 
through the normalization condition: 



S{t) + I{t) + R{t) = 1, 



(1) 



and they obey the following system of differential equa- 
tions 1311: 



dS_ 
dt 
cU 

It 
dR 



-q{k-l)IS, 

-III + q{k-l)IS, 



-j7 = 



(2) 
(3) 
(4) 



These equations can be interpreted as follows: infected 
nodes become recovered at a rate /i, while susceptible 
nodes become infected at a rate proportional to both the 
densities of infected and susceptible nodes. Here, q is the 
infection rate and k is the number of contacts per unit 
time. 

In the simulation, systems with size iV < 2^^ have been 
taken. Time is discretized and only the infection rate q 
is used as a parameter; infected people recover within 
one unit of time and can infect susceptible individuals 
with probability q connected directly to them within the 
same time scale. Initially all nodes are susceptible, one 
arbitrary node is chosen and taken to be infected. For 
the same network, 600 such choices have been taken and 
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FIG. 1: (Color online) Variation of Rsat with q for dif- 
ferent values of (5 = 0.0,0.5,0.9 and 1.0. Insets show the 
data collapse where Y = Raat has been plotted against 
X = {q-q,)N'/\ 
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FIG. 2: (Color online) Variation of Rsat with q for different 
values of (5 = 1.4 and 1.7. Insets show the data collapse where 
Y = RsatN'^''^ has been plotted against X = {q - qc)N^''' . 



quantities are averaged. A secondary averaging is made 
by considering 100 different network configurations. Dy- 
namics takes place in parallel; all the individuals revise 
their epidemiological state in one time step. Periodic 
boundary condition has been used in the simulation. 



III. RESULTS 
A. Static results 

1. Fraction of recovered population 

We calculate the fraction R{t) of recovered nodes as a 
function of time. This means the total fraction of nodes 
who are recovered till the time of measurement. It is also 
identical to the fraction of the population who have been 
infected at some point of time in the past. R{t) reaches a 
saturation value Rsat < 1 which depends on both q and 
5. We use synchronous dynamics to update the state of 
the nodes. The total duration r of the epidemic is also 
estimated. 

As a function of the infection probability q, Rsat is 
plotted for different values of 5 (Fig. [T]). From the knowl- 
edge that SIR shows the same criticality as the isotropic 
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percolation problem, it can be expected that a continuous 
phase transition takes place here. Since it is a finite size 
system, one can use the conventional finite size scaling 
method to estimate the critical (threshold) value qc and 
the associated critical exponents. Rsat vanishes below qc 
in the thermodynamic limit and should have a scaling 
form similar to magnetization. Therefore, the following 
finite size scaling form for Rsat is used: 



RsatO,N-^/'g,iiq-qc)N'^n. 



(5) 



where gi is a scaling function. We obtain data collapse 
with appropriate values of /3, qc and (Fig. [l]insets). Let 
us first discuss the region < S < 1 where qc — 0.430 
does not show any appreciable dependence on S. For S — 
0, the results should be comparable to a (addition type) 
Watts- Strogatz type network. The percolation threshold 
is known to be 0.401 there Our result for qc is 

slightly higher which could be because the average degree 
is different and also the fact that nearest neighbor links 
are always present here such that in comparison, long 
range neighbors are smaller in number. 

Results for < (5 < 1 show that the value of /3 ~ 
1 and ~ 3 are also unchanged in this region. This 
is not surprising, this region is known to have a small 
world behavior. Hence, mean field behavior is expected 
to be valid here and we find the value exponent P — I 
indeed matches with the mean field value for percolation. 
Assuming v = vd, where d = 6 is the upper critical 
dimension 0, @, IS, [l^l , one gets v = 1/2 from the fact 
that ~ 3, which also coincides with the mean field value 
for percolation phenomenon. 

For (5 > 1, (7c increases and approaches 1 at a value 
of 5 close to 2. This is expected, as for 5 > 2.0, the 
network behaves as a regular network in one dimension 
where qc — 1- The value of the exponent v however, 
appears to be unchanged while /3 shows a variation with 
(5, although not very strong. j3 starts deviating from the 
mean field value at around 5 = 1.4. 

The total duration r also shows a dependence on g; its 
peak value increases with the system size (Fig. |3]). This 
data is also analyzed by finite size scaling assuming the 
scaling form 



^N^"'g2{{q^qc)N^'% 



(6) 



where (/2 is another scaling function, u turns out to be 
very close to 3.0 from the above analysis as well for all 5 
values, while n shows a dependence on 5 ior 5 > 1.3. 

We show the values of the threshold values q^ and the 
exponents and j3 as functions of 5 in Fig. [5] Here 
qc shows variation with 5 for 5 > 1.0, the values of the 
exponents start differing from the mean field value at 
5 > 1.3. qc appears to reach 1 slightly above S — 2.0, we 
have checked that for 5 = 2.2, (jc is very close to 1. 

If we put q = qc, the variation of r with N shows that it 
diverges algebraically with N sublinearly in the entire re- 
gion < 5 < 2. This is interesting: r may be regarded as 
the minimum number of steps connecting two individuals 
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FIG. 3: (Color online) Variation of r with q for differ- 
ent values of (5 = 0.0,0.5,0.9 and 1.0. Insets show the 
data collapse where Z = tN~^^'^ has been plotted against 
X = (g - Q,)iVi/^ 
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FIG. 4: (Color online) Variation of r with q for different values 
of 5 = 1.4 and 1.7. Insets show the data collapse where Z = 
T-jy-f/"^ has been plotted against X = {q - qc)N'^^" . 



for exactly g = 1 in the network, i.e., it is comparable to 
the diameter of the network. However, one might expect 
that for q > qc, when spanning occurs, a finite number 
of the nodes can be reached by the infection procedure 
and hence r gives an estimate of the average number of 
steps connecting two individuals for the entire network 
(which may not necessarily be the minimum path). We 
know that in the small world region, the diameter scales 
as log(A^) but here we get a different scaling for r. We 
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FIG. 5: (Color online) The exponents and the transition point 
shown against S. 
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FIG. 6: (Color online) Frequency distributions for outbreak 
size Raat are shown for A'' — 2048 for two values of S. Data 
for q < qc,q > qc and q = qc are shown. For q < qc the 
distribution is unimodal while for g > gc, it is bimodal with 
two peaks occurring at Rsat ~ 1/A'' and Rsat ~ 0{1) (see 
text). Insets show that the data at g = gc has a power law 
variation. 



will discuss this point in section IV and V. 



2. Distribution of the outbreak size 

In the discussion in the preceding subsection, results 
for the average value of the outbreak size Rsat has been 
reported. We have also studied the distribution ^{Rgat) 
of Rsat which shows that it is bimodal in nature for q > qc 
and unimodal for g < i^c as found in ^34, 35]. ^(Rsat) 
has a peak near for all values of q which appears for 
the cases when the initially infected single node cannot 
transmit the disease to anyone else. The peak value at 
decreases with q as expected. For q > qc the $(i?sat) 
has a secondary peak at a larger value of Rsat ^ 0{1) 
(Fig. IHl). It may be noted that the peak values at Rsat ^ 
and Rsat ~ 0{1) are comparable for q > qc which 
makes the average value of Rsat a meaningful quantity. 
At qc, ^{Rsat) varies continuously and the tail of ^{Rsat) 
has power law decay as shown in the inset of Fig [51 The 
associated exponent varies with 5, for example, it is ~ 
1.41 for S = 0.5 and - 1.13 for 5 = 1.5. 

Such a bimodal behavior of distribution of outbreak 
sizes, which occurs as a single agent is assumed to be 
infected initially, is consistent with empirical data of dis- 
ease spreading analyzed by Watts et al [36|]. 
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FIG. 7: (Color online) The infected fraction of nodes shown 
against time for 5 = 1.5 for a system of size A'' = 2048. There 
is a peak at tp for q > qc. Inset shows the variation of total 
recovered fraction of nodes with time for 5 — 1.5. 



B. Dynamical results 

The growing density of recovered population R{t) has 
been plotted in the inset of Fig. [71 The data shows 
the expected saturation and a very fast growth at initial 
times suggesting an exponential behavior at early times 
[37| . It is found that the numerical data for R{t) can be 
fitted to the form: 



Rit) 



a exp(t/T) 

i+ccxp{t/T) ~ r 



(7) 



where a, c and T depends on q and S. The boundary 
condition assumed in the fitting is i?(0) = 0. 

The fraction of nodes infected at time t, I{t) = R{t + 
1)— i?(i) is plotted in Fig. [T] An initial growth and a peak 
value occurs at time t = tp only for q > qc- For q < q^ 
one gets a decaying behaviour right from t — 0. Such a 
decaying behaviour can occur if secondary infections are 
less than primary infections. Hence approximately. 
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which gives 



q<l/ik-l) 



(8) 



(9) 



Hence this argument can explain the absence of the peak 
for q < qc- The fact that the recovered population is no 
longer susceptible has been ignored in this argument, but 
for initial times, this will not matter when the recovered 
population is very small. This is supported by the data 
presented in Fig [S] which show that for q < qc, Rsat has 
very small values only. 

The variation of tp against q shows tp increases sharply 
from zero close to the transition point qc (Fig. [5]). Thus 
one can get an independent estimate of qc from this study. 
One may also plot the peak value I{tp) as a function of 
q which again shows an increase from zero close to the 
transition point (Fig. 

One can check how good is the fitting (eq. [7]) by the 
following approach. Assuming continuous time, I{t) can 
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FIG. 8: (Color online) The time at which the peak value 
occurs against q shows a sharp rise near qa (Data shown for 
N = 2048). 




FIG. 9: (Color online) The peak value of I against q shows a 
monotonic increase (Data shown for A'' = 2048). 



be expressed as Let Fit) = then from eq. ([7]), 

f exp(Vr) 



(10) 



(l + cexp(t/T))2- 

At large times, F[t) should fall exponentially with I{t) 
indeed shows such a behavior. The value oil = t^ where 
F(t) is maximum is found by solving the equation 

dF _ ^ eMt/r) - f&(exp(VT))2 _ 

dt (l + cexp(t/T))3 ' ^ ' 

which gives 



in 



Tin 



At tp = Tln(i), the value of F{t) is 



F{tp) 



dR 
'dt 



a 



(12) 



(13) 



Putting the values of a, c and T obtained from the fitting, 
tp and F{tp) can be computed from eqs p2)) and ^3]) 
and compared to the actual data. It shows very good 
agreement showing the quality of the fit. The comparison 
is shown in Table I for a particular value of q. 



TABLE I: a, c and T for three values of 5 and comparison of 
tp and F{tp) obtained from fitting and data(g = 0.58) 
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FIG. 10: (Color online) Variation of Rsat with q for mean 
field SIR model. 



IV. MEAN FIELD RESULTS FOR FINITE SIZES 

One can formulate a mean field type recursion relation 
with discrete time steps for the SIR model in which the 
average degree is (k) from eq. 

Iit)^qi{k)-l)Iit-l)S{t-l). (14) 

As I{t) = R{t+l)-R{t) and 5= l-R-I, this gives, 

R{t + 1) - R{t) = q{{k) - l){R{t) - R{t - 1))(1 - R{t)). 

(15) 

The initial conditions are i?(0) — 1/N and R{1) — 
R{0) + g(fc)i?(0)(l — R{0))- Obviously, here one assumes 
that the neighbors to which the infection spreads can be 
anywhere. We note that the system size enters through 
the initial conditions only, and it has been assumed that 
only one individual is infected in the beginning. 

R{t) and r are numerically estimated for different sizes 
using (fc) = 4 and shown in Figs [TU] and [11] The threshold 
value of q is very close to the theoretical estimate qc = 
l/((fc)-l) = l/3[l|. 

These data however do not show any finite size scaling 
behavior, i.e., a data collapse cannot be obtained. This 
is compatible with the fact that for a mean field model, 
finite size scaling does not work [s^. For the Euclidean 
model considered in this paper on the other hand, it is 
possible to obtain mean field exponents by finite size scal- 
ing in a finite region of the parameter space. 
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FIG. 11: (Color online) Variation of r with q for mean field 
SIR model. Inset shows the behavior of r-max with system size 
for mean field SIR model. 

Similarly, r also does not show finite size scaling. But 
exactly at q^, it has a peak value wrhich scales with N as 
log(iV). This agrees with the result of [3§|. Interestingly, 
such log(A^) behavior for r is not obtained aX q = qc 
for any value of 5 for the Euclidean network. This is 
consistent with the earlier results [13, EH that for small 
world networks, when percolation is considered, the giant 
component that emerges does not have a small world 
geometry. 



V. DISCUSSIONS OF THE RESULTS 

Before discussing the results, it is useful to recapitu- 
late briefly the studies already made on this network as 
comparison with these earlier works is necessary. 



A. Earlier studies on this network: Static 
properties 

It has been mentioned earlier that there are some re- 
sults available for the static properties of the network. 
A number of studies have been made to investigate the 
behavior of the network by calculating numerically the 
average shortest path as a function of the network size 
[i^ - lisj . Two earlier works [3 in which the greedy al- 
gorithm was used to evaluate the shortest paths, showed 
contradictory results. The study made with much larger 
networks indicated that the network had finite dimen- 
sional behavior for 1 < 5 < 2 [i^ while in [44], it was 
claimed that up to (5 = 2, the small world behavior ex- 
ists. In a more recent work '^ol], all possible shortest 
paths were evaluated numerically using a burning algo- 
rithm, and it was again found that the network retains 
the small world property vcp to 5 — 2 while the clustering 
coefficient vanishes below 5—1. Such a result was also 
obtained in |46j. However, one can still argue that it is 
an effect of finite sizes |45| . Simple scaling arguments 



suggest that there r nay be two transitions occurring at 

5 = I.O and 5 = 2.0 [3]. 

B. Critical phenomena studied earlier on this 
network 

The Ising model has been studied on the Euclidean 
network and again in two different studies contradictory 
results were obtained. While in [s?] it was found that 
the network behaves as a finite dimensional system for 
1 < (5 < 2, the study with the larger network [l^] in- 
dicated a small world behavior up to (5 = 2.0. Both 
studies, however, showed that the transition temperature 
remained 5 independent up to (5 = 1. 

SIR is expected to be related to percolation phenomena 
as already mentioned. Percolation has been studied in 
this type of network recently ^1] . The detailed study was 
made in two dimensions. In one dimension, results for 
two values oi 5 — 1.5 and 1.75 showed that the exponents 
are 5 dependent. The values of the threshold infection 
probability turn out to be larger than our results. This is 
probably because the networks are constructed in slightly 
different ways and also could be because site percolation 
was considered in [48i] . 

C. Present results and comparisons 

In comparison, the results obtained in the present work 
do not apparently lead to any definite conclusions for the 
behavior of the network for 1 < (5 < 2, while the mean 
field behavior obtained for 5 < 1 is consistent with the 
earlier results. In the mean field region, one has to use 
the upper critical dimension in the scaling relation to 
extract the exponents just like in the case of Ising model 

To analyse the region 1 < 5 < 2 carefully, the results 
for the static exponents v and /3 may be considered by 
comparing to the exponent values of isotropic percola- 
tion in different dimensions, v remains a constant in this 
region while (3 shows a change. It may be possible that 
v remains constant with both u and d varying. However, 
for percolation, we do not have such invariance of vd for 
dimensions 1-6 [i^. This suggests that the dimensional- 
ity is actually not changing. /? on the other hand does 
show a deviation from its mean field value (equal to 1) 
for 5 > 1.4. However, the deviation is not appreciable; 
the extrapolated value oi (3 aX 5 ~ 2.0 does not reach 
anywhere close to zero as it should. Rather, the value of 
/? for 5 even very close to 2.0 corresponds to what one 
would expect for a five dimensional lattice which indi- 
cates that in all probability it has a constant value up to 

6 — 2.{) but shows small changes in finite systems due to 
the transition occurring ai 5 = 2.0. 

One can also consider the dynamic exponent ji which 
shows a stronger dependence on 5 in the region 1 < 5 < 2 
compared to j3. Here it may be recalled that even in the 
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small world region, /i has a finite indicating t has a power 
law behavour while in the finite mean field case, r shows 
a logarithmic behaviour (sec IV) . Hence in contrast to the 
static exponents, /i does not reflect the behaviour of the 
network even in the small world region. Thus conclusions 
about the network's properties cannot be made on the 
basis of the behavior of ^. 

We have to compare the present results with those of 
static percolation on the same network which is available 
in [il]. For the results of the two points reported in this 
paper, the criticality appears to change appreciably as 
5 varies between 1 and 2 and does not show mean field 
behavior. The disagreement may be because of several 
factors: the construction of the two networks and the av- 
erage degree is different, secondly, it could be that the 
correspondence between percolation and SIR is not en- 
tirely true for the Euclidean network in this region. Such 
lack of correspondence has been noted in some earlier 
works Also, there could be such strong finite size 
effects that trying to compare results from two different 
studies is not useful. 



VI. CONCLUDING REMARKS 

From the discussion of the last section, we arrive at the 
following conclusion; SIR on this network shows mean 
field behavior up to (5 = 2.0 for at least the system sizes 
considered here. Evidently therefore, the network does 
not behave as a finite dimensional lattice for the region 
1 < (5 < 2 as has been claimed in some earlier works but 
effect of finite sizes are quite strong in the region 5 > 1.4. 

Talking about dimensionality, we would like to add a 
few words here. It is customary to use the upper critical 
dimension d and express i> — vd to get the mean field 
critical exponents using finite size scaling in these types 
of systems [2| . Hence one uses d = 4 for the Ising model 



and d = 6 for the percolation case. Surely one cannot 
assign an unique dimensionality to the network in this 
way at least in the mean field regime. However, assigning 
a dimension to the network in the mean field region is 
perhaps meaningless. 

What happens when mean field results are no longer 
valid? Theoretically, it is possible that one extracts a 
dimensionality in between 4 and 6 (which is what we are 
getting in the present work if we consider the value of 
/3) where percolation will not show mean field results but 
Ising model will. This appears far fetched and compels 
us to believe that the entire region < S < 2 has mean 
field behavior in agreement with [47[. 

The critical point shows variation with d in the in- 
termediate region as did the transition temperature for 
the Ising model [s^, 113 ■ only that, the time dura- 
tion also shows appreciable change in the scaling with 
the system size at criticality here. Apparently, the net- 
work having appreciable clustering leads to these changes 
which is felt in the intermediate region but overall small 
world property is retained such that mean field behavior 
prevails. 

We move over to a more profound question, do we reach 
any conclusion about what is the actual behavior of the 
network for 1 < 6 < 2? Still this issue is not resolved. 
As far as Ising model criticality and the SIR process is 
concerned, the latest results suggest it has a small world 
behavior. However, results of static percolation and some 
geometric properties indicate otherwise. The only possi- 
ble solution to justify these earlier results may be that 
at finite sizes, the network shows a small world behavior 
and deviations are only perceptible at very large length 
scales. 
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